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Abstract 

We write exact renormalization-group recursion relations for near- 
est-neighbor ferromagnetic Ising models on Migdal-Kadanoff hierar- 
chical lattices with a distribution of aperiodic exchange interactions 
according to a class of substitutional sequences. For small geometric 
fluctuations, the critical behavior is unchanged with respect to the uni- 
form case. For large fluctuations, as in the case of the Rudin-Shapiro 
sequence, the uniform fixed point in the parameter space cannot be 
reached from any physical initial conditions. We derive a criterion to 
check the relevance of the geometric fluctuations. 

PACS numbers: 05.50; 64.60; 75.10.H 
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There are a number of recent investigations on the critical behavior of 
ferromagnetic spin systems with the inclusion of aperiodic exchange interac- 
tions In particular, there are detailed studies of the critical behavior in 
the ground state of the quantum Ising chain in a transverse field 0, |3], [|, [5| 
(which is known to be related to the transition at finite temperatures of the 
two-dimensional Ising model). The nearest-neighbor ferromagnetic exchange 
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interactions are chosen according to some substitution sequences, and the ge- 
ometric fluctuations are gauged by a wandering exponent u> associated with 
the eigenvalues of the substitution matrix of the sequence. According to 
a heuristic criterion proposed by Luck[^], the critical behavior remains un- 
changed (that is, of Onsager type) for bounded fluctuations (small values of 
uj), but large fluctuations should induce much weaker singularities, similar 
to the case of a disordered Ising ferromagnet. 

In a very recent publication 0, we took advantage of the simplifications 
brought about by a Migdal-Kadanoff hierarchical (MKH) lattice |7j to per- 
form some exact calculations for characterizing the critical behavior of ape- 
riodic ferromagnetic Ising models. The (layered) exchange interactions be- 
tween nearest neighbors were chosen according to a certain class of gener- 
alized Fibonacci sequences @. In this paper, we review some of these cal- 
culations, discuss the well known case of the Rudin-Shapiro sequence, and 
present an exact derivation of an analog of Luck's criterion to check whether 
the geometric fluctuations are strong enough to change the critical behavior 
of the uniform system. 

Consider a particular two-letter generalized Fibonacci sequence given by 
the substitutions 

A^AB, 

B -> AA. 1 > 

If we start with letter A, the successive application of this inflation rule 
produces the sequences 

A -> AB -> ABAA -> ABAAABAB (2) 

At each stage of this construction, the numbers and Nb, of letters A and 
B, can be obtained from the recursion relations 



N' A 
N' b 




(3) 



with the substitution matrix 

m=i ; - ). in 




The eigenvalues of this matrix, Ai = 2 and A2 = — 1, govern most of the 
geometrical properties. For any sequence, the total number of letters, at 
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a large order n of the construction, depends asymptotically on A". The 
fluctuations are of order | A2 1 - Then, it is interesting to define the wandering 
exponent || , 

hi|A 2 | , . 

" = W (5) 

that expresses the asymptotic dependence of the fluctuations with the total 
number of letters, AiV (n) ~ N u . 

The nearest-neighbor Ising model is given by the Hamiltonian 

n = -^2J iJ a i a j , (6) 

(ij) 

with the spin variables a* — ±1 on the sites of a hierarchical diamond struc- 
ture. In Fig. 1, which is suitable for the period-doubling Fibonacci rule of 
Eq. ([!]), we draw the first stages of the construction of a diamond lattice with 
a basic polygon of four bonds (that is, of a MKH lattice with cell length 
b = 2, and number of branches q = 2, which amounts to 4 bonds in the 
diamond unit cell). As indicated in this figure, we simulate a layered system 
by the introduction of the interactions J a > and Jb > along the branches 
of the structure. From the rules of Eq. (|I|), it is straightforward to establish 
the recursion relations 

/ 2xa x b 
Xa ~ 1 + xW ' [ ) 



and 



x > = JE* ( 8 ) 



where xa = tanni^, i B = tanhKg, Ka = /3Ja, K b = /3Jb, and (3 is the 
inverse of the temperature. 

It should be remarked that similar procedures can be used to consider 
much more general substitutional sequences. However, to avoid any changes 
in the topology of the hierarchical lattice, we restrict the analysis to period- 
multiplying substitutions. In these cases, the largest eigenvalue of the in- 
flation matrix, Ai, gives the multiplication factor of the period. Therefore, 
b — Ai for all cases under consideration (and, in particular, b — Ai = 2, for 
the diamond lattice of the figure). 
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Fig. 1- Some stages of the construction of a Midgal-Kadanoff hierarchical lattice with 
bond length b — 2 and q = 2 branches (which is called a diamond lattice) for the period- 
doubling sequence A — > AB and B — > (letters A and B indicate the exchange interac- 
tions, J a and Jb). 



In the uniform case, J a 
recursion relation 



J 



B 



X 



J, Eqs. ( 

2x 2 
" 1 + a; 4 ' 



and (H) reduce to the simple 

(9) 



with two trivial and stable fixed points, x* = and x* = 1, and a non-trivial 
and unstable fixed point, x* = 0.543689..., which come from the polynomial 
equation 



x 



2aT + x 



x 



x 



2x + l 



0. 



(10) 



In the aperiodic case (J a ^ Jb) under consideration, the A components 
of the coordinates of the fixed points in the physical sectors of the xa — %b 
space (0 < Xa, x b < 1) come from the solutions of the equation 



Q 
A 



+ 2x\ - 4x 3 a + x A 



x a 



x A 



-2^4 + 1 



r 4 

x A 



+ 2 Xj 



1 



0. 



11' 



Therefore, a comparison with Eq. ( |T0| ) shows that the only fixed points are 
located along the xa = %b direction, and given by the same values as in the 



uniform 



x* B = 0, x* A = x* B = 1, and x* A = x* B = 0.543689.... The 



linearization about the non-trivial uniform fixed point yields the asymptotic 
equations 



Ax'a 
Ax\ 



CM 



B 




(12) 



where M T is the transpose of the substitution matrix, and the structure 
factor C is given by 

C = ~ Xa = 0.839286.... (13) 



x A 
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The diagonalization of this linear form gives the eigenvalues 



A 



i — 



CX 



i — 



2C = 1.678573..., 



(14) 



and 



A 2 = C\ 2 



C 



0.839286...., 



(15) 



where Ai = 2 and A 2 = —1 are the eigenvalues of the substitution matrix. 
As Ai > 1 and 1 7V2 1 < 1, the fixed point is a saddle node with a flipping 
approximation. Therefore, given the ratio t = J^j J a between the exchange 
interactions, the critical temperature is defined by the flow into this uni- 
form fixed point. From Eqs. (§) and fllPf), we see that the same eigenvalue 
Ai characterizes the (unstable) fixed point of the uniform model. Thus, in 
this particular example, with the wandering exponent u — 0, the geometric 
fluctuations are unable to change the critical behavior with respect to the 
uniform system. We can draw a phase diagram, in terms of the ratio r and 
the temperature T, where the critical line displays the same (universal) ex- 
ponents of the uniform case. Also, it is not difficult to check that the same 
sort of behavior (saddle point; largest eigenvalue associated with the uni- 
form system) still holds for all finite values of the branching number q of the 
diamond (b = 2) structure. 

To give an example with another value of the wandering exponent u, 
where the geometric fluctuations become relevant, consider the (four-letter) 
Rudin-Shapiro sequence[|], A — > AC, B — > DC, C -> AB, and D -»■ DB. 
The substitution matrix is given by 



Mrs 



I 1 1 \ 
11 
110 

V 1 1 / 



(16) 



with eigenvalues Ai = 2, A 2 = — A 3 = y/2, and A4 = 0, and the wandering 
exponent u — 1/2. For the Ising model on the diamond hierarchical lattice 
of the figure, it is easy to write the set of recursion relations 



Again, there are two trivial fixed points, x* A = x* B = x* c = x* D = and 1, and 
the non-trivial uniform fixed point, x* A = x* B = x* c = x* D = 0.543689..., as 
in the uniform case. The linearization about this uniform fixed point gives a 
matrix relation of the same form as Eq. (|i~2|) , 

(Ax') = CM T RS (Ax) , (19) 

with the same structure factor C, given by Eq. (|13|), and the eigenvalues 

Aj = CXx =2C = 1.678573..., (20) 

A 2i3 = CA 2 , 3 = ±CV2 = ±1.186930..., (21) 

and 

A 4 = CA 4 = 0. (22) 

Therefore, besides being unstable along the diagonal direction [xa = x B = 
%c = %d), this uniform fixed point is also unstable along two additional di- 
rections in the four- dimensional x A — x B — %c — x d parameter space. Given 
the ratios between the exchange interactions, there is no temperature associ- 
ated with any physical initial conditions in this parameter space so that we 
can reach the uniform fixed point. The critical behavior is of a much more 
subtle character as compared with the uniform case. 

Now it is interesting to devise an analog of Luck's criterion to gauge the 
influence of the geometric fluctuations. As we have seen in the previous 
examples, the largest eigenvalue in the neighborhood of the uniform fixed 
point is given by 

A 1 = AiC = bC, (23) 

where it is important to remark that the calculations are always performed 
for substitutional sequences and MKH lattices such that Ai = b. The second 
largest eigenvalue is given by 

A 2 = X 2 C = yA x . (24) 

Therefore, the fluctuations are relevant if 

|A 2 | = lM Al > 1. (25) 
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^From the exact recursion relations between the free energies associated with 
successive generations of a uniform ferromagnetic Ising model on a MKH 
lattice 0, we can write 

A 1 = b yt =b^ } (26) 



where a is the critical exponent of the specific heat of the uniform model 
and D is the fractal dimension of the lattice. From the definition of the 
wandering exponent, given by Eq. (|5|), we can also write 

|A 2 | = ff". (27) 

Inserting these expressions into inequality (0), we show that the geometric 
fluctuations are relevant for 

to > u c = 1 - — — . (28) 
I — a 

In the particular case of the diamond hierarchical lattice (6 = 2 and q — 2) 
the fractal dimension is given by 

so the criterion is reduced to the inequality 

a 

lo > u c = . (30) 

I — a 

As Ai = 2C = 1.678573..., we have a = -0.676533..., and to > 0.252764..., 
which explains the universal behavior of the first example {u — 0) and the 
relevance of the fluctuations in the case of the Rudin-Shapiro sequence (u> = 
1/2). This same criterion explains the change in the critical behavior of an 
aperiodic Potts model on the diamond hierarchical lattice above 4 + 2y / 2 
states, as recently shown by Magalhaes, Salinas, and Tsallis||. 
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Figure Caption 



Fig. 1- Some stages of the construction of a Migdal-Kadanoff hierarchical 
lattice with bond length b = 2 and q = 2 branches (which is called a diamond 
lattice) for the period-doubling sequence A — > AB and B — > AA (letters A 
and B indicate the exchange interactions, J a and J B ). 
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